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SUMMARY 


The service life of aircraft structural components undergoing random stress 
cycling was analyzed by the application of fracture mechanics. The initial crack 
sizes at the critical stress points for the fatique-crack growth analysis were estab- 
lished through proof load tests. The fatique-crack growth rates for random stress 
cycles were calculated using the half-cycle method. A new equation was developed 
for calculating the number of remaining flights for the structural components. The 
number of remaining flights predicted by the new equation is much lower than that 
predicted by the conventional equation. 


INTRODUCTION 


Aircraft structural components commonly contain flaws , defects , or anomalies of 
various shapes; these either are inherent in the basic material or are introduced 
during the manufacturing and assembly processes. A large percentage of service 
cracks found in aircraft structures are initiated from crack nucleation sites such 
as tool marks, manufacturing defects, and surface microinclusions (ref. 1). Under 
the combined influences of environment and service loading, these flaws may grow 
to reach catastrophic sizes, resulting in serious reduction of service life or com- 
plete loss of the aircraft. Thus, to a great extent, the integrity of the aircraft 
structure is dependent upon the safe and controlled growth of cracks as well as the 
achievement of residual strength in their presence. 

The operational life (service life) of aircraft structural components is 
affected by the magnitude and cumulative effects of external loads coupled with 
any detrimental environmental action. The presence of moisture, chemicals, sus- 
pended contaminants, and naturally occurring elements such as rain, dust, and sea- 
coast atmosphere can cause deterioration in structural strength due to premature 
cracking and acceleration of subcritical crack growth (refs. 2 and 3). 

As aircraft structures begin to age (that is, as flight hours accumulate), 
existing subcritical cracks or new cracks can grow in some high-stress points of 
the structural components. The usual approach is to inspect the structures peri- 
odically at certain intervals. However, even after inspection there may be some 
undetected cracks in a structure. To ensure that the structure still has integrity 
for future flights, proof load tests are usually conducted on the ground. The pur- 
pose of proof load tests is to load a structure to certain proof load conditions 
(slightly lower than the design limit load conditions) to test its integrity. If 
there should exist undetected cracks in the structural component that are larger 
than the critical crack sizes, that structural component will surely fail during 
the proof load tests and will be replaced. This process can reduce the chance of 
catastrophic accidents during flight. If all structural components survive the 
proof load tests, then fracture mechanics can be applied to estimate fatigue life 
(number of remaining flights) for each critical structural component by using the 
initial crack size established for each structural component during the proof load 
tests and then using the stress cycles (obtained through strain gage measurements) 
for each structural component during the first flight after the proof load tests. 


This report describes the application of fracture mechanics and the half- 
cycle method to calculate the number of remaining flights for aircraft structural 
components . 

The author gratefully acknowledges contributions by William W. Totton and 
Jules M. Ficke of Synernet Corp. in setting up computer programs for flight data 
reductions and crack growth computations. 


SYMBOLS 


A crack location parameter 

a crack length of edge crack, one-half the crack length of through- 

thickness crack, or depth of surface crack 

a^ crack size (length or depth) after the £th flight 

a° limit crack length (or depth) associated with the operational peak load 

P 

a c initial fictitious crack length (or depth) established by the proof 

load tests 

C material constant in Walker crack growth rate equation 

c half-length of surface crack 

E complete elliptic function of the second kind 

Fj£ number of remaining flights after the £th flight 

Fi number of remaining flights after the Ith flight predicted from the 

newly developed equation of remaining flights 

f operational peak stress factor (f < 1) 

i integer associated with half-cycles, or critical stress 

points 

j integer associated with half-cycles 

Kj Mode I stress intensity factor 

Kj c Mode I critical stress intensity factor 

K ma x Mode I stress intensity factor associated with o max , AMKO max "Vna/Q 
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k modulus of elliptic function 

l integer associated with flights 

M k flaw magnification factor 

m Walker exponent associated with K max 

N number of stress cycles available for operations 

N£ number of stress cycles used during the ith flight 

n Walker exponent associated with stress ratio R 

Q surface flaw shape and plasticity factor 

R stress ratio, a min /a max 

t thickness of plate 

V A front hook vertical load 

V BL left rear hook vertical load 

v BR right rear hook vertical load 

Aa^ amount of crack growth during the £th flight 

6ai crack growth increment resulting from one cycle of constant-amplitude 

stress cycling under loading magnitude of A and Ri 

Ak Mode I stress intensity amplitude, AM K (a max - a m i n ) V^a/Q 

remote uniaxial tensile stress 

00 

o° uniaxial tensile stress associated with the peak operational load level 

p 

proof tensile stress induced by the proof loads 
ck tensile stress at critical stress point i 

* 

peak value of Oi induced by the proof loads 
c max maximum stress of the stress cycle 

o m _^ n minimum stress of the stress cycle 
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tensile strength 


yield stress 

ultimate shear strength 

<J> angular coordinate for semielliptical surface crack 


THEORY 


Fracture Mechanics 

The top part of figure 1 shows the most common types of cracks: through- 

thickness crack, surface crack, and edge crack. According to fracture mechanics, 
the stress intensity factor Kj for the Mode I deformation (tension mode) associated 

with any type of crack can be expressed as 



(1 ) 


where A is the crack location parameter (A = 1 for the through- thickness crack, 

A = 1.12 for both the surface and the edge cracks (see fig. 1))? M^ is the flaw 
magnification factor (for a very shallow crack M^ - 1; as the depth of the crack 
reaches the back surface of the plate, Mk = 1*6 (see fig. 1)); is the remote 

uniaxial tensile stress; a is one-half the crack length for the through- thickness 
crack, or the length of the edge crack, or the depth of the surface crack (see 
fig. 1); and Q is the surface flaw shape and plasticity factor given by 


Q = [E(k)]2 



( 2 ) 


where is the yield stress and E(k) the complete elliptic function of the second 
kind, defined as 


E (k ) 



1 - sin^ <(> d<)> 


(3) 


where is the angular coordinate for a semielliptical surface crack, defined in 
figure 1 , and the modulus k of the elliptic function is defined by 

k -V 1 - 

where c is the half-length of the surface crack (see top center of fig. 1). 

The bottom part of figure 1 shows the plots of a/2c as functions of Q for dif- 
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ferent values of o^/o^ and the plot of M k as a function of a/t, where t is the 
plate thickness. 


Proof Load Tests 

The purpose of proof load tests is to load the entire aircraft structure (or 
its components) to certain proof load levels to test structural integrity and to 
establish initial fictitious crack sizes associated with critical structural com- 
ponents for fatigue life analysis. The proof load levels are usually slightly 
lower than the design limit load conditions associated with different maneuvers. 

If there exists in a certain structural component a previously undetected crack 
that is larger than the critical crack size associated with the proof load, that 
component will certainly fail during the proof load tests and will be replaced. 

Thus, a catastrophic accident during flight can be avoided. If the entire struc- 
ture survives the proof load tests, then the critical stress point of the struc- 

p 

tural component has been subjected to a proof tensile stress a induced by the 
proof loads. If Kj c denotes the critical stress intensity factor (or material 
fracture toughness) of the structural component material, the maximum crack length 

p 

a c the structural component can carry under the proof loads without failure (or 
rapid crack extension) may be calculated from equation (1) by setting Ki ■ Kj , 

P c 

and = o^. In reality, there may not be any cracks that developed during proof 

p 

tests; however, it is assumed that a fictitious crack of length a c has been cre- 
ated at the critical stress point of the structural component during the proof load 
tests. During actual operations, the structural component will be subjected to much 

P o 

lower stress levels than the proof stress level o m , if is defined as the peak 
operational stress level (highest peak of the stress cycles), then according to 

equation (1), the structural component can carry a fictitious crack of size a°, 

P o 

which Ls much larger than a c . The value a c thus determined is considered to be 

P 

the limit crack size toward which the initial crack a c is allowed to grow after 

o P 

repeated operations. The crack size difference, a c - a c , is then the crack size 

increase permitted for the structural component in repeated operations. The left- 

hand plot in figure 2 shows crack length, a, as a function of normalized stress 

o /a , where a is the tensile strength of the material. It is seen that the lower 
o° / u U 

the operational stress level, the larger the limit crack size available for the struc- 
tural component. The right-hand plot in figure 2 shows crack length, a, as a function 
of number of constant stress cycles N. 


Remaining Flights 

If the structural component is cycled under constant stress amplitude (an ideal- 
ized case for the purpose of discussion) for N-j cycles during the first flight with 
the associated crack growth of Aa-j , then the number of remaining flights F-j before 
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the limit crack size is reached may be estimated from the following conventional 
equation (see fig, 2): 


o P 
a c “ a c _ N 

Aa-j ~ Nt 


( 5 ) 


Equation (5) is based on the assumption that the amount of crack growth per flight 
for all subsequent flights will be equal to Aa-j of the first flight. In reality, 

the amount of crack growth per flight will steadily increase with the number of 
flights accumulated, and the actual number of remaining flights will be less than 
the value F-j predicted by equation (5) if the number of remaining flights is large 
(that is, F-j >> 1), As discussed in the following section, for a relatively low 
range of F-j , equation (5) may give a reasonably accurate prediction of the number 
of remaining flights. The amount of crack growth Aa-j in equation (5) may be calcu- 
lated from the following Walker equation (refs. 4 to 6) for fatigue-crack growth 
rate under constant-amplitude stress cyclings: 


= C ( K max ) m ( 1 - R) n = C(AKP(1 - R) n “ m (6) 

aN 

where C, m, and n are material constants, and K max , Ak, and R are respectively maxi- 
mum stress intensity factor, stress intensity amplitude, and stress ratio, given by 


K max 


AM K a max 



(7) 


Ak - ^K^max ” a min) 



( 8 ) 


R ~ a min/ a max 


(9) 


where o max and o m j, n are respectively the maximum and the minimum stresses of the 
constant-ampli tude stress cycles. 

In reality, the stress cycles encountered during operations at the criti- 
cal stress points of the structural component are not constant-amplitude stress 
cycles. To apply equations (6) to (9) to variable-amplitude stress cycles, dif- 
ferent methods must be developed. In this report a half-cycle method (refs. 7 to 
10) is used in the calculation of the fatigue-crack growth rate for variable- 
amplitude stress cycles. 


Half-Cycle Theory 

The top part of figure 3 shows an example of random stress cycles (variable- 
amplitude loading history). The stress history curve is the combination of a 
series of both increasing and decreasing load half-cycles of different loading 
magnitude (Ak, R) , as shown in the lower part of figure 4. The half-cycle theory 
(or half-wave theory) (refs. 7 to 10) states that the damage (or crack growth) 
caused by each half-cycle of either increasing or decreasing load is assumed to 
equal one-half the damage caused by a complete cycle of the same loading magnitude 
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(Ak, R). This means that the damage caused by the complete cycle could be equally 
divided between the two phases of increasing and decreasing loads. The loading 
sequence thus can be resolved into half-cycle groups of increasing and decreasing 
loads (see lower part of fig. 3). Each half-cycle (either increasing or decreasing 
load) can then be considered as a half-cycle of the constant-amplitude cyclings 
under the same loading magnitude (Ak, R) and can be computed separately in time 
sequence to estimate the corresponding damage. The half-cycle theory thus permits 
accurate evaluation of the load spectrum from a recorded load time history. If 

(£ « 1 r 2, 3, . . . ) is the final crack length after (£ = 1, 2, 3, % . .) random 

stress cyclings in the £th flight, then according to the half-cycle theory, a^ may 
be calculated from 


with 


2N 


^ 6a.; 


a £ = a £ 


i=1 


2N fl 


6a 


Aa £ = Z ~Y- - a* - a Jl _ 1 and a,,., = a, 


i=1 


( 10 ) 


(11 ) 


where 6a±/2 is the crack growth increment induced by the ith half-cycle under the 
loading magnitude Ak^ and ; 6a^/2 is assumed to equal the crack growth increment 

induced by a half-cycle of the constant- amplitude stress cycle fatigue test under 
the same loading magnitude (Ak^, R^). Thus, by using equations (6) to (9), 6a£/2 
may be calculated from 


6a-L 
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n r m n-m 

Rj.) * "2 ( AKi ) (1 - Ri) 


( 12 ) 


The crack length, a, associated with AK^ (see eq. (8)) will be the summation of the 

initial crack length and all the crack growth increments created by all the pre- 
vious ha If -cycles : 


r> 6a~; 

a - 4 + I (i » 2) 

j»1 


(13) 


Similar to the case of constant-amplitude stress cycling (see eq. (5)), the number 
of remaining flights (£ = 1, 2, 3, . • .) after the £th flight of random stress 

cycling may be calculated from 


F * = 


a c - a A-1 
Aa a 


a c “ a *_1 


a c ~ 
2N„ 


1 


i=1 


a £ " a Jl-1 


6a i 

2 


(14) 



Figure 4 graphically illustrates how to evaluate the crack growth increment Sa^/2 

associated with the ith half-cycle of the random stress cycling by using the plots 
of AK| as functions of da/dN for different values of R. 

As mentioned previously, equations (5) and (14) both may give relatively 
accurate values for F£ when the number of remaining flights is relatively low. 

When predicting large F^ , equations (5) and (14) must be modified because Aa£ 
increases with the flights accumulated. 


New Equation for Remaining Flights 

For the case of constant-amplitude stress cycling, the amount of crack growth 
Aa £ ( £ = 1 , 2, 3, . , . ) during the £th flight (N^ cyclings) may be obtained by 

intearatina equation (6): 


Aa, 


i - C^AM K a ma x 
where equation (7) was used. 


kin /, «\ n m/2 

(1 - R) M a *-i) 


( i = 1, 2, 3, . . .) 


(15) 


For simplicity, if a max , R» and remain the same for all the flights, then 
the following crack growth rates can be established by using equation (15), assuming 
*V a i-i << 1: 

m/2 


Aa : 
Aa i 


P \m/2 . 

a c + Aai \ m Aai 

r ] — 1 t • • • 

P ) 2 P 

a c / a c 


(16) 


Aai 

Aai 


Aa<; 

Aai 


a 3 


m/2 


m/2 


P . A \ m /2 

a c + Aai + ^ a 2 

T 

a c 


= 1+ 2 7 


P A . \m/2 

a c + Aa-| + Aa 2 + Aa 3 ' 

~P~ 
a c 


Aai 

\ 

P 

a c , 

r 


^m ^ a 1 

• 3 ( 

2 P 
V a c , 


(17) 


(18) 


Aai 


m/2 


£-1 


a' + Aai + ^ a 2 + ^ a 3 + 
a c 


• + Aa 


£-1 


\m/2 


= i + U-i )| 


m 

2 P 


(19) 
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If the available crack size, a c - a c , can allow F-j number of remaining flights, then 


F i terms 


1 ° - a ? Aa 1 + ^ a 2 + ^ a 3 + • • • + Aa^ + . . . + Aaj^ 


( 20 ) 


Aai 


Aai 


where the left-hand side is F-j , the number of remaining flights predicted by assuming 
that Aai = Aa 2 =* Aa 3 = . . . = Aa^ = . . . * Aap 1 (see eq. (14) for i - 1). Sub- 
stituting equations (14) and (16) to (19) into equation (20), 

F •) terms (Pi - l) terms 

' K N m Aai f A “ I N 

Fi=(1+1+1+. . .+1) + j — [1 + 2 + 3 + . . . + (Fi - 1 )] (21) 


or 


m Aa! (Fi - 1 ) [(Fi - 1 ) + 1 ] 
F 1 = Fi + 2 P 2 


( 22 ) 


Equation (22) may be rearranged into the form 


Aai 

— F 1 
a c 


( m Aa-i \ _ 

1 '“TT" r 1 

V 4a c / 


- Fi 


Solving for Fi , 



which gives the relationship between Fi and Fi . 


(23) 


(24) 


If the prediction of the number of remaining flights is based on the crack 
growth Aa^ that occurred during the £th flight, equation (24) takes the form 


iW /7E5 

ra Aa J t \V a *-i 


- 1 


(25) 


Equations (24) and (25) both apply to the case of constant-amplitude stress cycling. 

However, they may be used for the case of variable-amplitude stress cycling without 

introducing significant error. Figure 5 shows the plot of equation (24) (that is, 

- P 

F-j as a function of Fi ) for m = 3.6 (Inconel 718 alloy) and Aai/a c * 0.01814 for 
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the example described in the following sections. Some typical values of F 1 and F-j 
are compared in the following tabulation: 


F 1 •• 1 

10 

50 

100 

150 

200 

500 

1000 

F 1 : 1 

9 

33 

53 

70 

84 

147 

218 


The ratio F -j /F-j is 0,53 for F-j * 1 00 and decreases to 0,218 for F^ = 1000 , Thus, 

equation (5) certainly overpredicts the number of remaining flights , and safety 
factors ranging from 2 to 4 must be used depending on the range of the value of F*j , 


EXAMPLE PROBLEM 


For this report, the example problem chosen for the fatigue analysis using the 
half-cycle method is the severe fatigue problem of the three hooks of the NASA 
B-52-008 carrier aircraft pylon used to carry the space shuttle solid rocket booster 
drop test vehicle (SRB/DTV) shown in figure 6 . The 49,000-lb SRB/DTV was attached 
to the pylon with one front hook and two rear hooks (see fig, 6 ). The SRB/DTV was 
carried up to high altitude and released to test the performance of the solid rocket 
booster main parachute. The shapes of the front hook and the two rear hooks are 
shown in figures 7 and 8 , respectively. The front hook is made of Inconel 718 alloy 
and the two rear hooks of AMAX MP35N alloy (AMAX Specialty Corporation). Table 1 
shows material properties of the two alloys. Because of the great weight of the 
SRB/DTV, the three hooks had serious fatigue life problems. Fracture mechanics and 
the half-cycle method can be applied to predict the service life (number of remaining 
flights) of the three hooks that carried the SRB/DTV. Reference 10 presents the 
detailed fatigue analysis of the NASA B-52 aircraft pylon major components. 


Critical Stress Points 

Before conducting fatigue-crack growth analysis, the locations of the critical 
stress points for the three hooks had to be determined. This was done by performing 
NASTRAN finite-element stress analysis of the three hooks (ref. 11). The critical 
stress point of each hook is located at the inner circular boundary of the hook. 
Figures 7 and 8 show the exact locations of the critical stress points of the three 
hooks. Through the NASTRAN stress analysis, the relationships between the stress 
at the critical stress point and the hook loads were established as 


«1 

= 7.3522 

X 

10- 3 

V A 

for 

front hook 

(26) 

°2 

= 5.8442 

X 

10-3 

VBL 

for 

left rear hook 

(27) 

0 3 

= 5.8442 

X 

IQ " 3 

Vbr 

for 

right rear hook 

(28) 


where o^, o 2 , and 03 are respectively the stresses (in kips, or 10 3 lb, per square 
inch ( ksi ) ) at the critical stress points of the front hook, left rear hook, and 
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right rear hook, and V A , V BL , and V BR are the corresponding hook vertical loads in 
pounds. During proof load tests and during flight, V A , V BL , and V BR were measured 

by means of strain gages installed near the critical stress points of the hooks 
(refs. 10 and 11). Equations (26) to (28) were used to generate stress cycles for 
the fatigue-crack growth analysis using the strain-gage-measured values of V A , V BL , 
and V BR . 


Initial and Operational Crack Sizes 


During the proof load tests, the three hooks were loaded to their respective 

P 

peak proof loads to establish the initial crack size a c for each critical stress 

* 

point. The peak proof stresses (i = 1, 2, 3) at the critical stress points 
induced by the peak proof hook loads may be calculated from equations (26) to (28). 

P 

The proof c rack size a c at the critical point of each hook established by the proof 
load tests may be calculated from equation (1) by setting Kj = Kj c : 



(29) 


where A = 1.12 for the surface crack, M K = 1 (which was obtained from the lower 

right plot of figure 1 for a/t << 1 because the depth of the crack is very small 
compared with the depth of the hook), and the value of Q, the surface flow shape and 
plasticity factor, will be determined as follows. 


The surface crack is assumed to be semielliptical in shape with an aspect ratio 

of a/2c = 1/4. (This value is based on the observation of surface cracks of the 

★ 

fractured old rear hooks.) Taking a/2c = 1/4 and o^/ay = a i/ a Y = 1 (because the 
growth of plastic zones around the critical stress points was neglected, the values 

of a i calculated for the three hooks slightly exceeded the corresponding yield 
stresses Oy), the curve for a^/Oy = 1 in the lower left plots of figure 1 gives 
Q = 1-25. 


* 

If the peak value of the stress cycles during operation (or flight) is ftfj_, 
where f is the operational peak stress factor (f < 1), then the operational limit 
crack size may be calculated from 


o 

a c “ 


Q 

TT 


K Ic \ 2 

AM K fcTi/ 



f2 


(30) 


* 

Part of table 2 shows the peak proof hook loads, the peak proof stresses and 

p 

the proof crack sizes a c at the critical stess points. Note that for all three 

★ 

hooks, a exceeded the failure stresses ay of the hook materials (see table 1), yet 

the three hooks did not fail during the proof load tests. The reason is that equa- 
tions (25) to (27), established by the NASTRAN analysis, are for a purely elastic 
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case without consideration of plastic deformations. In reality the plastic zone 
can develop around the critical stress point , and therefore the hook can actually 
carry a greater load than the brittle failure load. In the present fatigue-crack 
growth analysis, only the elastic case is considered. 


Load Spectra 

To perform the fatigue-crack growth analysis, the load spectra (stress cycles) 
for the three critical points of the hooks must be obtained first. Using the 
strain-gage-measured values of the three hook loads, V A , V BL , and V BR during 

the first test flight, the three stresses o$_ (i = 1, 2, 3) may be calculated by 

using equations (26) to (28). Figures 9 to 1 1 (ref. 10) show portions of the 
loading histories (load spectra) calculated for the critical stress points of the 
three hooks during a takeoff run. Those load spectra were obtained by filtering 
the original load spectra down to 5 Hz to eliminate the small-amplitude high- 
frequency stess cycles that are considered unimportant in the present fatigue 
life analysis. Notice that the load spectra for all three critical stress 
points exhibit a high degree of random cycling. 


Calculations of Crack Growth 

To apply the half-cycle method, the load spectra (see figs. 9 to 1 1 ) were first 
resolved into a series of half-cycles of increasing and decreasing loads of differ- 
ent loading magnitude (Ak, R) (see fig. 3). The crack growth increments Sajy2 per 

half -cycle were calculated from equation (12) in time sequence and summed (using 
eq. (10)) to give the total amount of crack growth per flight for each critical 
stress point. Finally, equations (14) and (25) were used to calculate the number 
of remaining flights associated with each critical stress point. It must be empha- 
sized that in using the half-cycle method, every half-cycle of different stress 
amplitude in the load spectrum is calculated, and thus the half-cycle method can 
give an accurate evaluation of the fatigue-crack growth as compared with, say, the 
exceedence-count method (for which some of the stress peaks lying below the mean 
line could be missed). 


Results 

Figures 12 to 14, taken from reference 10, show the fatigue-crack growth curves 
calculated for the three critical stress points for the first test flight; the 
maneuver transition points are indicated. Note that for the three hooks, the 
fatigue-crack growth rate is greatest during the initial stage of taxiing and 
the takeoff run and becomes very low during cruising because of relatively low- 
amplitude stress cyclings. Table 2 lists the amount of crack growth Aa-j , opera- 
tional peak stress factor f, operational crack size a°, and remaining flights F^ 
and F-j calculated respectively from equations (14) and (24). Figure 15 shows the 

plots of F<| and Fi as functions of f for the three hooks. Note that the over- 
prediction of remaining flights based on Fi (eq. (14)) becomes more pronounced as 
the number of remaining flights increases. The arrows in figure 15 indicate the 
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actual operational peak load levels. At these load levels, safety factors must be 
in the range of 2 to 2.5 if F-j is used instead of F-j . 


CONCLUSION 


Fracture mechanics and the half-cycle method were applied to the service life 
analysis of aircraft structural components. The initial crack sizes at the crit- 
ical stress points of the structural components were determined by using proof load 
tests. The random stress cycle fatigue-crack growth rates were calculated using 
the half-cycle method. A new equation was developed for calculating the number of 
remaining flights for the structural components. The newly developed equation pre- 
dicted the number of remaining flights more accurately (a much lower number) than 
did the conventional equation (which is based on the assumption that the amount of 
crack growth per flight remains constant). 
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Figure 1 . Three types of cracks and the plots of surface flaw shape 
factor and flaw magnif ication factor. 
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-Equation (24) 


m * 3.6 (Inconel 718) 
Aa« /a? ■ 0.01814 


Figure 5 . Number of remaining flights predicted from Ko equation 
of remaining flights . 
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Figure 6 . Geometry of space shuttle solid rocket booster 
drop test vehicle (SRB/DTV) attached to B~52 pylon (view 
looking inboard from right side). 
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Figure 13. Crack growth curve for stress point 2 ( 02 ), flight 1. 



Figure 14. Crack growth curve for stress point 3 ( 03 ), flight 1 
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